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Abstract 

Epitaxial CaMn03 films grown with 2.3% tensile strain on (OOl)-oriented LaA103 substrates are found 
to be incipiently ferroelectric below 25 K. Optical second harmonic generation (SHG) was used for the 
detection of the incipient polarization. The SHG analysis reveals that CaMn03 crystallites with in-plane 
orientation of the orthorhombic b axis contribute to an electric polarization oriented along the orthorhombic 
a (resp. c) axis in agreement with the predictions from density functional calculations. 

PACS numbers: 77.55.Nv, 77.80.bn, 71.15.Mb, 42.65.Ky 
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I. INTRODUCTION: MULTIFERROICS WITH STRAIN-DRIVEN FERROELECTRICITY 



Control of the magnetic response by electric fields and of the dielectric response by magnetic 
fields is highly interesting for spintronics applications and devices based on a rigid coupling of 
magnetic and dielectric properties. The possibly most fertile source for magnetoelectric cross 
correlations are compounds with a coexistence of magnetic and electric long-range order, called 
multiferroics. 1 ' 2 The magnetoelectric phase coexistence is inherently rare: The established mech- 
anisms driving magnetic and ferroelectric order require partially filled and empty 3d orbitals, 
respectively, and are therefore mutually exclusive unless they arise from different lattice sites 
This lead to an intense search for alternative mechanisms promoting ferroelectricity in the pres- 
ence of long-range magnetic order. A variety of mechanisms were identified up to now. 5 One of 
the most notable examples is BiFe03, the only robust ambient multiferroic, in which ferroelectric- 
ity is driven by an electronic 6s lone pair of the Bi 3+ ions. Furthermore, in compounds where the 
magnetic long-range order breaks the inversion symmetry and induces a spontaneous polarization 
the magnetoelectric effects are intrinsically strong because of the rigid coupling of the improper 
ferroelectric polarization to the proper magnetic order parameter(s).- 1 ^"- 

In spite of the impressive magnetoelectric coupling effects observed in some of these com- 
pounds an ambient multiferroic with pronounced and strongly coupled magnetization and polar- 
ization has not been discovered thus far. As a consequence, the ambition to discover fundamentally 
new mechanisms promoting multiferroic order is unbroken. In the past, new multiferroics were 
usually obtained via the chemical route by synthesizing materials from chemical building blocks 
with a high potential to promote multiferroic order. As alternative it was recently suggested to 
take a structural route and modify the unit cell parameters of non-multiferroic compounds until 
they develop the magnetoelectric phase coexistence. An efficient way to achieve this is strain. The 
improvement in thin-film growth technologies in the past years allows the growth of thin films 
with epitaxial strains of several percent.— Such strain can be used for increasing ferroelectric-^ 
and ferromagnetic^ 2 - 1 ^ transition temperatures, induce ferroelectricity,— or enhance the magneti- 
zation of a ferromagnet.— In 2010 it was shown that tensile biaxial strain of 1.1% pushes EuTi03 
into a multiferroic state by inducing a spontaneous polarization estimated as 29 /iC/cm 2 that co- 
exists with a ferromagnetic phase. — First-principles density functional calculations explained this 
transition by softening of a polar phonon mode driven by spin-phonon coupling.— 

Aside from this landmark experiment, strain-induced ferroelectricity complementing magnetic 
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order was also predicted for EuO and CaMn03 . For EuO, epitaxial tensile or compressive strain in 
the order of as much as 4% is required 18 which is difficult to achieve experimentally. CaMn03 is 
a particularly attractive candidate compound. First, the threshold tensile epitaxial strain required 
for stabilizing a ferroelectric phase was predicted to be a moderate 2.1%.— Second, multiferroic 
CaMn03 would create an interesting exception to the aforementioned 3d rule because in contrast 
to other multiferroic perovskite oxides antiferromagnetism and ferroelectricity in CaMn03 would 
both be associated to the Mn 4+ cation site. At the bulk level, the antiferrodistortive structural 
instability is calculated to dominate the ferroelectric instability, so that the ground-state is non- 
polar. However, epitaxial tensile strain is expected to promote the polar mode to the extent that 
the compound can eventually become spontaneously polarized in thin films. Third, "chemical 
strain" exerted by replacing Ca by Sro.55Bao.45 in bulk crystals already revealed the emergence of 
a ferroelectric state.— 

Unfortunately, the experimental verification of ferroelectricity in strained CaMn03 films re- 
mained an unsolved task thus far. For CaMn03 a polarization of 4 /iC/cm 2 was predicted to 
emerge (at K) along the orthorhombic a axis for tensile strain of 2.1% applied in the orthorhom- 
bic ac plane.— This value is orders of magnitude larger than the spontaneous polarization ob- 
tained in the magnetically induced ferroelectrics. 7,8 However, pyroelectric current measurements 
are spoiled by leakage currents and by the inefficiency of generating electric poling fields and 
pyro-currents within the film plane. A more suitable experimental approach for the verification of 
in-plane ferroelectricity is clearly required. 

Here we report the presence of an incipient ferroelectric state below 25 K in epitaxial CaMn03 
films subjected to 2.3% tensile strain. Incipient ferroelectrics are compounds with emerging 
electric-dipolar long-range order where the Curie-Weiss temperature 7cw is so close to zero (either 
Tew ?t or 7cw ^ 0) that the long-range-ordered state is not yet stabilized down to K.— "— In- 
cipient ferroelectrics with Tew ^ are sometimes termed "quantum paraelectrics"— "— when the 
suppression of the ferroelectricity is related to quantum fluctuations and distinguished from the 
"incipient ferroelectrics" in a narrower sense. However, here we follow the majority of published 
work and refrain from making this distinction. One of the best known incipient ferroelectrics is 
SrTi03 (Tqw > O)^ 1 ^ so that we use it as reference compound to which we compare our CaMn03 
data. 

The polarizability of the CaMn03 films is demonstrated by optical SHG. The direction of the 
incipient polarization is derived from symmetry considerations and found to be oriented diagonally 
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between the principal in-plane axes associated with the cubic perovskite subcell of the pseudocubi- 
cally grown films. This orientation agrees with earlier predictions derived from density functional 
calculations; here a refined density functional approach backs up our conclusions. 



H. PROBING FERROELECTRICITY BY SHG 

SHG is a well established tool for probing ferroic order in bulk crystals and thin films.— ^ The 
nonlinear optical process describes the generation of a light wave at the frequency 2 ft) in a material 
with ft) as the frequency of the incident light.— This is described by the expression Pi(2co) = 
EoXi jkE j(oo)E k (co) . The component Xijk of the nonlinear susceptibility tensor couples j and k 
polarized contributions of the electric field E( ft)) of the incident light to an i polarized contribution 
to the polarization P(2ft>) driving the SHG light field. In the electric-dipole approximation % is 
a polar tensor so that components Xijk 7^ are obtained in non-centrosymmetric systems only.— 
Thus, SHG is well suited for detecting ferroelectric order breaking the inversion symmetry.— In 
contrast to linear optical techniques the ferroelectric SHG signal emerges free of background at 
the Curie temperature. SHG is particularly useful for probing leaky or in-plane ferroelectricity 
because in contrast to pyro-current measurements the finite conduction does not interfere with the 
detection of the spontaneous polarization and electrodes are not applied. In addition, the degree 
of freedom of spatial resolution inherent to optical techniques allows one to probe the spatial 
distribution of the spontaneous polarization and, hence, image domains. 



HI. PEROVSKITE SUBCELLS IN PSEUDOCUBIC SAMPLES 

Bulk CaMn03 crystallizes in a distorted orthorhombic structure described by the space group 
Prima. .— In orthorhombic coordinates the lattice parameters are a = 5.279 A, b = 7.448 A, and 
c = 5.264 A. This structure possesses a perovskite subcell that can be approximated with a/y/2 ~ 
b/2 « c/y/2 ~ ^cube — 3.72 A as cubic lattice parameter as shown in FigJTJ Below — 122 K 
CaMn03 exhibits a G-type antiferromagnetic order. According to Ref. 19, a polar ground state 
is obtained for tensile strains larger than 2.1% for which a cubic substrate with £? C ube 

> 3.80 A is 

required. 

Our CaMn03 films (thickness 40 nm) were grown on (001) -oriented LaA103 substrates with 
a LAO = 3.81 A as lattice constant of the pseudocubic subcell. This corresponds to 2.3% tensile 
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strain of the films. Pulsed laser deposition was used for the epitaxial growth of the CaMn03 films. 
The substrates were kept at a constant temperature of 650 °C during the deposition, which was 
carried out at a pressure of 0.04 mbar of flowing oxygen. After the deposition, the samples were 
cooled to 400 °C maintaining the same conditions. The oxygen pressure was then increased to 
300 mbar, followed by slow cooling to room temperature. The structural study was carried out by 
x-ray diffraction (XRD) using a Seifert XRD 3000P for the 0-20 scans (Cu K a , A = 1.5406 A). 
The films were shown to be homogeneous and the structure corresponds to the composition of the 
target (Ca/Mn = 1) in the limit of accuracy. Sharp and intense diffraction peaks (see Fig. Old)) 
suggest neatly crystallized single-phase films. Using the XRD results, the out-of-plane parameter 
of the films was calculated to be < 3.71 A, confirming that the films are under epitaxial tensile 
strain. 

Figure |2] shows a transmission electron microscopy (TEM) image of one of the CaMn03 films. 
The TEM image shows that the film is composed of nested regions with a lateral size of < 10 nm. 
Fourier transformation (FT) reveals that there are three different types of regions that are exemplar- 
ily highlighted by colored circles. The corresponding states are associated to the three orientations 
of the orthorhombic unit cell of CaMn03 with respect to the substrate, i.e., with the orthorhombic 
b axis pointing along the x, y, or z direction of the substrate lattice. Applying this assumption, 
those spots in the FT data that are uniform across the sample are related to a simple cubic per- 
ovskite cell. The remaining spots can be explained by applying the distorted CaMn03 perovskite 
structure with the same rotations of the Mn06 octahedra that are present in the bulk. As detailed 
in Fig.[2tb) the corresponding three sets of spots are then identified as (h k/2 0), (h/2 k 0) and (h/2 
k/2 0). They are associated to [101] (sets 1 and 2) and [010] (set 3) zone axes patterns of the bulk 
orthorhombic Prima structure with cell parameters a = \[2a c , b = 2a c , c = \/2a c . Thus, on the one 
hand, the substrate enforces its cubic lattice parameters onto the CaMnC>3 film, but on the other 
hand, the CaMn03 retains the bulk orthorhombic atomic distortion for each of its three orientation 
states. 

We recall, that within the pseudocubic approximation the orthorhombic a axis and c axis may 
be interchanged so that we have a total of six possible orientations for the orthorhombic CaMn03 
cell on the LaA103 substrate. However, since the resolution of the TEM experiment does not allow 
us to distinguish between the a axis and the c axis, only three different orientations are identified 
in Fig. [2] 

Note that the TEM study indicates the presence of a small amount of secondary phase con- 
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eluded from the observation of supplementary spots along the (110) direction of the perovskite 
subcell but not detectable in the XRD data. Such superstructure is compatible with the reduced 
phases of CaMnOs — and is possibly related to the strained state of the film. 

IV. SHG ON FERROELECTRIC CaMn0 3 

Knowing the appropriate framework for the description of our CaMn03 films, i.e. the pseu- 
docubic approximation, we can now derive the possible symmetries of the ferroelectric phase and 
the resulting polarization selection rules for SHG. 

Non-polar CaMn03 is centra symmetric with mmm as orthorhombic point symmetry. In the 
pseudocubic approximation we neglect the difference between the a and the c axis which changes 
the point symmetry to 4 /mmm with the b axis as the fourfold axis. The spontaneous polarization of 
the strained pseudocubic unit cell may be oriented parallel (case i) or perpendicular to the b axis. In 
the latter case it may be oriented along the principal a or c axis (case ii) or diagonally, including an 
angle of 45° with these axes (case iii). Cases (i) and (ii, iii) reduce the point symmetry to 4mm and 
mm2, respectively, with '4' and '2' indicating the direction of the spontaneous polarization. Lower 
symmetries do not have to be considered because they would correspond to other, unphysically 
arbitrary directions of the spontaneous polarization. 

In Table U the SHG contributions for cases (i) to (iii) are given for all the possible orientations 
of the orthorhombic unit cell within the pseudocubic lattice. Only tensor components Xijk that can 
be addressed with light incident perpendicular to the CaMn03 film (k\\z) are considered. This 
excludes all the components with i,j, or k = z since this would involve longitudinally polarized 
light. Note that the net SHG yield obtained from the CaMn03 film is a mixture of the SHG 
contributions for all the orientations of the orthorhombic unit cell that are possible within the 
pseudocubic lattice. 

V. EXPERIMENTAL SETUP 

Prior to our SHG experiments, we performed pyroelectric current measurements on the 
CaMn03 films. For this purpose, two gold electrodes with a gap of about 1 mm were grown 
onto the surface of the films. However, as in previous experiments the leakiness of the epitaxial 
CaMn03 films in combination with the in-plane geometry spoiled the polarization measurement. 
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This leaves SHG as approach for the detection of a spontaneous polarization. 

In the SHG experiments, frequency-tunable laser pulses of about 130 fs are emitted from an op- 
tical parametric amplifier which is operated at 1 kHz by a Ti: sapphire amplifier system. The SHG 
data are taken in the spectral range 2hco = 1.8 — 3.0 eV which covers the lowest O 2 — > Mn 4+ 
charge-transfer excitation and the Mn 4+ (?2 g ) — > Mn 4+ (e g ) intraband transfer.—"— In order to sup- 
press any surface-induced SHG contributions not coupling to ferroelectric order a near-normal- 
incidence reflection geometry with a reflection angle incident to the surface normal of approxi- 
mately 2° is employed. The polarization of the incident laser beam is set by a half-wave plate. 
Optical filters are chosen to suppress background light in front of the sample and to separate the 
SHG wave from the fundamental light behind the sample. The polarization of the SHG light is ana- 
lyzed with a Glan-Taylor prism. For further spectral filtering the signal light is transmitted through 
a monochromator. It is then detected by a photomultiplier tube and normalized to a reference SHG 
signal in order to account for spectral variations of the fundamental light and of the efficiency of 
the SHG setup. In addition, the spectra were normalized by dividing the SHG signal obtained 
from the sample by the spectrally flat reference SHG response of a silver mirror. Alternatively, 
for obtaining spatial resolution, the samples were imaged onto a liquid-nitrogen-cooled camera 
chip using a standard telephoto lens with a resolution of about 25 /im. A liquid-helium-operated 
cryostat is used to cool the samples to temperatures between 4.5 and 300 K. Electric fields were 
applied to the sample via polished steel plates with a diameter of about 1 cm in between which the 
samples were mounted. 

VI. EXPERIMENTAL RESULTS AND DISCUSSION 

Prior to the experiments on the CaMn03 films the SHG response of CaMn03 and LaA103 sin- 
gle crystals was investigated in order to identify any bulk background contributions. As mentioned 
before LaA103 is a centra symmetric insulator with a band gap of 5.6 eV— 9 and a pseudocubic per- 
ovskite subcell that exhibits no ferroic order, so that electric-dipole-type SHG contributions are 
not expected. In agreement with this, a SHG signal from pure LaAlC>3 was not observed in the 
temperature range from 5 to 300 K and for photon energies between 1.8 and 3.0 eV. A CaMn03 sin- 
gle crystal grown by the floating-zone method revealed a spectrally flat, temperature-independent 
SHG signal that was recorded between 4.5 and 150 K. It may be related to defects, surface contri- 
butions, or SHG contributions beyond the electric-dipole approximation. 
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Figure |3] shows the normalized spectral, temperature, and polarization dependence of the SHG 
signal obtained on the epitaxial CaMn03 films. The pronounced spectral dependence of the SHG 
signal at 5 K demonstrates that it is not related to the background contribution observed on the 
CaMn03 bulk crystals. The spectral dependence is independent of the tensor component. A peak 
of the SHG intensity is observed at around 2.4 eV, which is probably related to the aforementioned 
hg — > eg transition of the Mn 4+ ion.— — For the following experiments a SHG energy of 2.1 eV 
was chosen because of the high intensity of the fundamental laser beam in combination with a 
reasonably large SHG yield. 

Figure [3ja) shows the temperature dependence of the SHG signal at 2.1 eV. Above 25 K a con- 
stant SHG background is obtained that shows an isotropic polarization dependence and is present 
up to at least 150 K. This is the background signal already observed on the bulk sample and can 
therefore be regarded as "zero bias" of the SHG measurement. At 25 K a pronounced SHG signal 
emerges and increases continuously towards 5 K. As Fig.Hd) reveals, the SHG signal of CaMn03 
has a distinct polarization dependence. The anisotropy measurement shows the SHG yield polar- 
ized parallel to the polarization of the incident fundamental light while rotating this polarization 
by 360°. This leads to four equally long lobes with maximum SHG intensity along the (110) 
in-plane-diagonal directions of the pseudocubic lattice. A fit entered as solid line into Fig. |3{d) 
reveals that the polarization dependence of the SHG signal is perfectly described by a single SHG 
tensor component: Xppp = Xaaa ( see Table U). This uniquely points to CaMn03 crystallites with 
an out-of-plane orientation of the orthorhombic b axis and a spontaneous polarization along the 
orthorhombic a (resp. c) axis, see Fig.djc). With two possible in-plane orientations of this po- 
larization the SHG signal from crystallites of either orientation adds up to reveal a fourfold SHG 
anisotropy in spite of the mm! symmetry of the individual crystallites. This conclusion is in per- 
fect agreement with the predictions made by Bhattacharjee et ah— and in Section IVII Al In both 
cases, density functional theory (DFT) is used to consider an out-of-plane orientation of the or- 
thorhombic b axis, and a spontaneous in-plane polarization along the orthorhombic a (resp. c) axis 
is found. 

Hence, both the SHG data in Fig. [3] and the DFT results point to an electric polarization in the 
pseudocubic epitaxial CaMn03 films that is induced by tensile lattice strain of 2.3%. In contrast to 
the bulk CaMn03 ground state, the polar mode becomes unstable in the strained film so that a polar 
displacement can appear. Note that the emergence of the SHG signal below 25 K cannot be due 
to the magnetic order. First, magnetization-induced SHG also be observed on the CaMn03 bulk 
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sample, and its polarization would have to be different from the polarization of the SHG signal in 
Fig. (2d). Second, DFT predicted that the antiferromagnetic order and its critical temperature of 
Tn = 122 K are not affected by the substrate strain and the resulting ferroelectric transition.— We 
thus conclude that strained pseudocubic CaMnC>3 constitutes a strain-driven multiferroic below 
25 K. 

In order to investigate the rigidity of the polar state, its response to electric poling fields and 
thermal cycling was investigated with the results shown in Figs.@]and[5l Figure |4] shows the effect 
of a static electric field on the polarization dependence of the SHG signal. The anisotropy of the 
SHG signal in the absence of an electric field and in a static field of 10 6 V/m during and/or after 
cooling from 40 to 5 K is compared. The electric field was applied along as well as diagonally in 
between the direction of the pseudocubic axes. (Only the latter case is depicted here since both 
cases lead to the same result.) According to Fig. |4] the electric field has no effect on the SHG 
yield. In addition, cycling the electric field between ±10 6 V/m showed no sign of polarization 
reversal or a hysteresis. 

Applying an electric field is expected to reverse the polarization of the domains with a polar- 
ization component antiparallel to the field and, thus, reduce the number of oppositely polarized 
domains, ideally towards a single-domain state. Since the SHG signal from oppositely polarized 
domains interferes destructively,— the field poling is therefore expected to enhance the SHG yield. 
Yet, the insensitivity of the SHG signal to the electric field shows that polarization switching does 
not occur. This may indicate that the polarized regions are so strongly pinned by the substrate 
strain and the pseudocubic twinning that the applied field is too low for polarization reversal. Yet, 
it is unlikely that even in a field-cooling experiment no tendency at all for the alignment of the 
spontaneous polarization along the applied field is observed, in particular in view of the large 
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value of the spontaneous polarization (and the related field energy) expected from DFT (Ref. 
and Section IVnA]). 

The spatially resolved SHG intensity of the sample is shown in Fig. \5\ The images reveal a 
grainy distribution of the SHG intensity in the form of resolution-limited speckles. According to 
Fig. Ob) the position and relative brightness of the speckles does not change when a consecutive 
annealing cycle through 60 K is applied. (Note that the arrangement of the speckles in the exem- 
plary red circles is the same in Figs. |5Ja) and|5Ib).) The result does not change when an electric 
field is applied during the temperature cycle. 

Figure|5]thus reveals that the insensitivity to electric-field cycling is matched by the insensitivity 
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to thermal cycling. Grainy distribution of the SHG intensity are a result of the interference of 
SHG contributions from areas with a size below the optical resolution limit. Here, these areas can 
either correspond to the differently oriented crystallographic regions constituting the pseudocubic 
structure in Fig. |2]or to a distribution of nanometer- sized ferroelectric domains. The similarity 
of Figs. [5ta) and [5^b) points towards the former. Domains would change in the course of an 
annealing cycle unless they are strongly pinned. However, as argued before, pinning effects that 
are pronounced enough to withstand an electric poling field when crossing the Curie temperature 
in the course of a temperature cycle are unlikely. 

We therefore conclude that the polarization picked up by the SHG signal is related to incipient 
ferroelectricity in the strained CaMn03 films. As we will see in the following this conclusion is 
supported by both experiment and theory, (i) In Fig. |3jc) we show the temperature dependence 
of the SHG signal from SrTi03 that we measured for comparison. The retarded emergence of the 
signal in this incipient ferroelectric is qualitatively very similar to the temperature dependence of 
the SHG signal from CaMn03. It is therefore obvious to expect incipient ferroelectricity to be also 
present in CaMn03. (ii) The SHG signal in SrTiC>3 displays a distinct spectral and polarization 
dependence (not shown) just like the SHG signal in CaMn03. (iii) As confirmed by Fig. |3jc) 
incipient ferroelectric s can display a pronounced SHG signal although there is no spontaneous 
polarization breaking the inversion symmetry. (We assume that with the high polarizability char- 
acteristic for incipient ferroelectric s the fundamental light wave itself drives the AC polarization 
that is breaking the inversion symmetry and probed by SHG. This was also proposed for explain- 
ing forbidden Raman lines in KTa03. 41 ) In the absence of magnetic or structural transitions in our 
CaMn03 films at 25 K, incipient ferroelectricity is thus a compelling reason for the emergence of 
the SHG signal, (iv) The proximity to the critical strain predicted by DFT— is in conformity with 
a ferroelectric potential that is too shallow to stabilize a spontaneous polarization. 



VII. VERIFICATION BY DENSITY FUNCTIONAL THEORY 



As mentioned before the observation of an incipient state with an in-plane polarization along the 
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However, in Ref. 
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orthorhombic a or c axis agrees well with the DFT results reported in Ref 
investigations were limited to CaMn03 films with an out-of-plane orientation of the orthorhombic 
b axis as shown in Fig. EDc). Since regions with an in-plane orientation of b are present in Fig. |2] 
we have to expand our DFT analysis now. We do this in two ways: (i) by verifying the results in 
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Ref. 19 using an alternative approach for their derivation; (ii) by adding the scenario of an out-of- 
plane orientation of b. In both cases we retain the pseudocubic approximation of the orthorhombic 
unit cell. 



A. Out-of-plane orientation of the orthorhombic b axis 

The formerly used pseudopotential approach with the generalized-gradient-approximation 
(GGA) Wu-Cohen exchange correlation functional (as it is implemented in the ABINIT pack- 
age) was replaced by the projector- augmented wave (PAW) methodology and the GGA PBEsol^ 
exchange correlation functional (as it is implemented in the VASP code.-^ii) Converged results 
were achieved with a plane-wave cutoff of 500 eV and a fc-point grid of 4 x 2 x 4. All calcula- 
tions were performed with collinear magnetism by taking the G-type antiferromagnetic order as 
the magnetic ground state. Phonon frequencies have been calculated for Pnma space group and 
different values of epitaxial strain using the frozen-phonon method with small atomic displace- 
ments of ±0.03 A. The cell parameter was imposed to be cubic and fixed to that of the substrate 
in two directions while it was relaxed in the third direction. 

First, we calculated the evolution of the soft transverse-optical (TO) modes versus the epitaxial 
strain applied along the a (resp. c) direction (see Fig.[TJc)). The strain dependence of the square 
of the frequency Q. of the three soft TO modes (TO a , TO^, TO,;.) polarized along the three direc- 
tions of the crystal is shown in Fig. [6] Interestingly, all three modes become softer under tensile 
epitaxial strain. Within the PBEsol functional the TO c mode becomes unstable (Q. 2 < 0) at a criti- 
cal epitaxial strain of 3.7%. This is followed by the TO a mode which becomes unstable at 4.0%, 
whereas the TO^ mode remains stable up to 5.0% , but with a low frequency. This clearly repro- 
duces the fact that CaMn03 develops a ferroelectric instability under tensile epitaxial strain that 
is resulting in an in-plane polarization. However, we note that the present calculations predict the 
occurrence of the ferroelectric instability at a critical tensile epitaxial strain of 3.7%. This value is 
larger than the value of 2.0% derived previously— and refers to the c instead of to the a axis. This 
difference can be related to the different functionals (GGA PBEsol instead of GGA Wu-Cohen) 
and the different DFT methodology (PAW instead of pseudopotentials) employed here. 

In order to scrutinize the influence of the choice of functional, we also performed calculations 
by the local density approximation (LDA) based on the PBEsol volume: We froze the volume 
and the cell parameters to the values obtained within the PBEsol functional and performed LDA 
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atomic relaxations and phonon calculations. In Fig. [6] the according values of £l 2 are shown for 
the three soft TO modes. We observe the same epitaxial strain sensitivity as reported with the 
GGA-PBEsol functional. However, now the ferroelectric instabilities develop at a lower epitaxial 
strain: 3.2% for the TO c mode, 3.6% for the TO a mode, and 4.8% for the TO/, mode. The result 
illustrates the considerable influence the choice of the functional has on the the critical epitaxial 
strain obtained for the emergence of ferroelectricity. Nevertheless, all calculations show the same 
trend, i.e., a ferroelectric a- or c-axis instability dominates the response of the CaMn03 lattice 
under tensile epitaxial strain applied along the a and c directions. 

The ground state was determined by condensing the unstable modes in the structure once they 
emerge. We then performed an atomic relaxation for each mode and checked the presence of re- 
maining instabilities. We thus found that up to a tensile strain of 5.0% the ground state is always 
obtained by the condensation of the TO c mode alone within the PBEsol and the LDA function- 
al. Hence, the condensation of the TO c mode removes all other instabilities. The ferroelectric 
polarization is then predicted to occur along the in-plane c direction. The amplitude of the polar- 
ization depends on the value of the epitaxial strain. With the Berry phase method 45 we obtain a 
spontaneous polarization of 12 /iC/cm 2 with the PBEsol functional and 18 /iC/cm 2 with the LDA 
functional for an epitaxial strain of 4.0%. 
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to identify individual atomic contribu- 



Performing the same analysis as in Table II of Ref. 
tions to the eigendisplacements of the unstable ferroelectric mode, we also find that the unstable 
TO c mode is dominated by Mn and O atomic motions (23% for Mn, 76% for O and only 1% for 
Ca at an epitaxial strain of 4% with the PBEsol functional), thus supporting the idea that ferroelec- 
tricity in CaMn03 is mainly driven by the Mn atoms at the B-site, in spite of their partially filled 
J-orbitals.— 



B. In-plane orientation of the orthorhombic b axis 

Applying the strain along the orthorhombic b direction (see Fig. CEJb)) constrains b to the v-cell 
parameter of the substrate. The a and c directions are not parallel to the substrate which constrains 
their projection to the the x-cell parameter of the substrate. However, the angle between a and 
c can relax since the out-of-plane projections of a and c are not constrained. 4 ^ For catching this 
scenario, we performed the calculations on a cell with 40 atoms. Figure [7] shows the evolution 
of Q. 2 for the three soft TO modes denoted TO x (mode polarized along the x direction of the 
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substrate), TQ^ (mode polarized along the y direction of the substrate with a small component 
along z), and TO vz (mode polarized along the y and the z directions of the substrate. The TO modes 
become unstable at critical epitaxial strains of, respectively, 3.6%, 3.8% and 3.9% within the LDA 
functional. With the PBEsol functional only the TO x mode becomes unstable at 4.0% whereas the 
TO-yM mode and the TO >z mode are still stable at this strain. In analogy to the behavior of the TO c 
mode with epitaxial strain along a and c, condensing the TO x mode for strain along b stabilizes 
the TO^ mode so that the ground state is given by the condensation of the TO A mode alone. At 
an epitaxial strain of 4.0% this gives rise to a spontaneous polarization of 4.6 /iC/cm 2 along the 
in-plane x direction with the LDA functional and a spontaneous polarization of < 0.01 /iC/cm 2 
with the PBEsol functional. 

While the critical epitaxial strain for inducing a ferroelectric instability depends on the 
exchange-correlation functional and the methodology used, all calculations confirm the SHG re- 
sults: Independent of the scenario (out-of-plane orientation of b with strain along a and c or 
in-plane orientation of b with strain along b), a spontaneous polarization along the orthorhombic 
b direction of the crystal is never stable. Furthermore, a polarization along the x-direction with 
in-plane orientation of b is not observed experimentally. This is consistent with the tendency for 
the higher values of the epitaxial strain that are required for its occurrence. 

Performing the same analysis as for the case where the strain is applied along the a and c 
directions, the eigendisplacement of the unstable TO A mode is also found to be strongly dominated 
by the Mn and O motions: 45% for Mn, 54% for O and only 1% for Ca at an epitaxial strain of 
4% (PBEsol functional). 

vm. CONCLUSION 

In summary, epitaxial pseudocubically twinned CaMn03 films with 2.3% tensile strain were 
found to be incipiently ferroelectric below 25 K. The polarizability was detected by optical SHG. 
According to the symmetry analysis of the SHG signal CaMn03 crystallites with out-of-plane 
orientation of the orthorhombic b axis contribute to an emerging polarization directed along the 
orthorhombic a (resp. c) axis in agreement with a variety of DFT approaches. The antiferromag- 
netic order as well as the incipient ferroelectric order arise from the Mn 4+ cation site which thus 
creates a remarkable exception to the 3d rule for perovskite multiferroics. With our experiments 
we emphasize the potential of strained perovskite oxides as resource for a rich variety of multi- 
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ferroic compounds. We expect that further research will either lead to constituents with a lower 
threshold for strain-induced ferroelectricity or to higher values of substrate-induced strain stabi- 
lizing the ferroelectric state beyond the incipient behavior found in the present CaMn03 films. 
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TABLE I: SHG contributions for all the possible orientations of the orthorhombic b-axis and the resulting 
spontaneous polarization P sp within the pseudocubic lattice. Case (i): P sp of the strained unit cell is oriented 
parallel to the fourfold fr-axis. Case (ii): P sp is perpendicular to the ft-axis and along the principal a or c 
axis. Case (iii): P sp is perpendicular to the &-axis, including an angle of 45° with a and c (these directions 
are denoted as o and p with o ± p). Only tensor components Xijk tnat can be addressed with light incident 
perpendicular to the CaMn03 film (k\\z) are considered. This excludes all the components with or k = z 
since this would involve longitudinally polarized light. 
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FIG. 1: Orthorhombic unit cell of CaMn03 and possible orientations of CaMn03 grown on LaA103. (a) 
Orthorhombic unit cell of CaMn03. (b, c) Principal orientations of CaMn03 grown on LaA103 LaA103 
with (b) in-plane and (c) out-of-plane orientation of the orthorhombic b axis. The orthorhombic axes are 
denoted as a, b, c whereas the pseudocubic axes are denoted as x, y, z, with a^o as the pseudocubic unit 
cell parameter, (d) 6-26 XRD scan of a typical strained CaMn03 film epitaxially grown on LaA103. The 
peak at 44° is caused by the sample holder. 
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FIG. 2: (a) High resolution TEM image of a CaMn03 film (thickness 40 nm) revealing regions with a differ- 
ent crystallographic orientation and a lateral extension of a few nanometer. On average, the image possesses 
a fluctuating contrast suggesting the presence of strain fields due to the nested configuration of the differ- 
ently oriented regions.—^ (b) Fourier transformation obtained from three different areas, corresponding 
to the three types of orientations of the orthorhombic CaMnOs unit cell within the film. The schematic 
drawing summarizes the four different sets of reflections that can be obtained in the Fourier transform of 
the perovskite subcell (white) and the three orientation variants (color). In total six orientations of the or- 
thorhombic CaMnOs unit cell are possible, since a and c axis may be exchanged within the pseudocubic 
approximation. 
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FIG. 3: (a) Temperature dependence of the SHG signal at 2.1 eV. At 25 K a pronounced polarized SHG 
signal emerges and increases continuously towards 5 K. Above 25 K we find a temperature independent, 
spectrally featureless SHG background that is present up to at least 150 K. (b) SHG spectrum of Xzzz at 5 K. 
The resonance at 2.4 eV is probably related to the t 2g — > e g transition of the Mn 4+ ion which is also observed 
in x-ray absorption spectra of bulk CaMn03.— (c) Temperature dependence of the SHG signal obtained 
from incipiently ferroelectric bulk SrTi03 for comparison. The similarity to the temperature dependence of 
SHG from CaMn03 is striking, (d) Polarization dependence of the SHG signal at 5 K and 40 K. The fit 
entered as a solid line is perfectly described by a single SHG tensor component: Xppp = Xaaa (see TableU). 
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FIG. 4: Electric-field dependence of the SHG anisotropy. Samples were zero-field cooled (ZFC, E = 0) or 
field cooled (FC, E = 10 6 V/m) from 40 K and measured at 5 K and 2. 1 eV with or without the electric field 
applied. The insensitivity of the SHG signal to the electric field indicates that no polarization switching 
occurs. Lines are fits to the data according to Fig. 0c). 
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FIG. 5: Spatially resolved SHG intensity of a CaMn0 3 film at 5 K and 2.07 eV SHG photon energy. A 
grainy distribution of SHG intensity in the form of resolution-limited speckles is obtained. The relative 
position and brightness of the speckles does not change after two consecutive cooling cycles through 60 K. 
The red circles mark exemplary areas. 
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Tensile strain (%) 

FIG. 6: Calculated square of the TO frequency Q. of CaMn03 versus a-c epitaxial strain. TO f( , TO;,, and 
TO f are the soft TO modes polarized along the orthorhombic a, b, and c axis, respectively. At a critical 
epitaxial strain of 3.2% (LDA functional) the TO f mode becomes unstable and initiates a net polarization 
along the orthorhombic c axis (pseudocubic diagonal) reaching 18 /^C/cm 2 at an epitaxial strain of 4%. 
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FIG. 7: Calculated square of the TO frequency CI of CaMnOs versus b epitaxial strain. The x, y and z 
directions are the directions of the polarization of the modes according to Fig.[TJb). The TO Y mode becomes 
unstable at a critical epitaxial strain of 3.6% (LDA functional) resulting in a net polarization of 4.6 /iC/cm 2 
along the x direction of the crystal (pseudocubic axis) at an epitaxial strain of 4%. 
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